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Recently it has been suggested that an enhancement in the visibility of ghost images 
obtained with thermal light can be achieved exploiting higher order correlations HI. This 
paper reports on the status of an higher order ghost imaging experiment carried on at 
INRIM labs exploiting a pseudo-thermal source and a CCD camera. 
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1. Introduction 

Quantum Imaging is a new quantum technology addressed to exploit properties of 
quantum optical states for overcoming limits of classical optics^. 

Various protocols have been recently proposed | 3 | 4 | 5 | 6 | 7 | 8 ] and realized 
| 9 | io | ii | i2 | i3 | i4 | ran g m g f rom super-resolution to sub-shot noise imaging. 

One of the oldest ideas, that can find very interesting applications is the so 
called Ghost Imaging. In summary, the idea at the basis of this protocol is that 
whether one disposes of two noise-correlated light beams, one crossing an object to 
be imaged and then detected by a bucket detector without any spatial resolution, 
the other addressed to a spatial resolving detector (as a CCD camera), then the 
image of the object can be reconstructed by considering correlations between the 
two measurements. 

This protocol was firstly predicted ^ and then demonstrated with twin beams 

and later with thermal light l 18 l 19 l 20 l 2 3 (with only a smaller visibility) . 

Recently, it was shown that the use of higher order correlation functions can 
improve the visibility of ghost imaging E 

A first experiment in this sense was realized in^l This interesting result, nev- 
ertheless, was valid only in the approximation of very high light flux where G^ 
was reconstructed by a G^ measurement and was realized with a " fake" mask (i.e. 
by artificially blackening some pixels of the camera via software). Very recently a 
second experiment appeared where a real object (a double slit) was imaged through 
a scan that allowed to reconstruct G^ 
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Fig. 1. Third order imaging scheme. There are two reference arms and a test arm. Distances d\ 
and d,2 have to satisfy the imaging condition. The detectors on the reference arms are spatially 
resolving while in the test arm a bucket detector is present. 



Here wc present our preliminary results toward third order correlation Ghost 
Imaging reconstruction of a real object observed through a CCD camera and com- 
pare the achieved results with second order ghost imaging ones. 

2. Theory 

The imaging scheme (Fig. [I]) counts three arms, two reference arms and a test 
arm. The unknown object is inserted in the latter. Each distinct imaging system 
is characterized by its impulse response function hi(x i ,Xi) } where i = 1,2,3, x i 
are the starting transverse coordinates and xi are the coordinates at the detection 
plane. Hence, considering a field Ei(x i ), its value at the detection plane is given by: 

Ei(xi) = J dx'jh^x'^x^E^x'i) (1) 

In our scheme the subscript i = 1 refers to the test arm. 

The detectors in the reference arms are spatially resolving detectors and record 
the intensity distribution J 2 (x2) and ij^xs) in the pixel x% and x% respectively. 
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Information about the object is obtained by measuring the (normalized) second 
order correlation coefficient c2 in the case of second order ghost imaging and the 
(normalized) third order coefficient c3 in the case of third order ghost imaging, 
defined as: 

c2 = < hjx^IjjXj) >-<I 1 (x 1 )>< Ij[Xj) > ^ 

\/% (zi)-\/M2(z.j) 

with j — 2,3 and: 

„/ v < (h(xi)- < h(x!) >)(h(x 2 )- < h{x 2 ) >){I 3 {x 3 )- < I 3 (x 3 ) >) > 

c3{x 1 ,x 2 ,x 3 ) = 3 r~T^ a^^^ 3^^^ 

VM3(»i) V/W2) vM^) 

(3) 

where < Ii(xi) >=< E* (xi)Ei(xi) > is the mean intensity of the i-th beam. In 
formula ^ we define: 

M2(»i) =< (Ii(xi)- < h{xi) >f > (4) 
as the second central moment, while in Q 

H3(Xi) =< (Ii{Xi)- < h{Xi) >) 3 > (5) 

is the third central moment. Symbols < .. > refers to temporal averages. As a 
matter of fact it is the correlation between the beams which is responsible of the 
image formation in the ghost imaging process EEH, for this reason we considered the 
(normalized) correlation coefficient c2 and c3 instead of the G^ and G^ functions. 
In fact, this is helpful as it removes unwanted effects related to the unbalancing of 
the intensity or of the fluctuations between the beams which may somehow hide 
the correlations of interest. 



3. Experimental set-up 

In our experiment (Fig. ([2])), the imaging system consists of three classical correlated 
thermal light beams. The thermal field is obtained by passing a coherent source 
through a random media. The source is a passively Q-switched microchip laser 
with a pulse duration of about 1 ns and a maximum average power of 1 W; this 
source together with a frequency doubler (the second harmonic generator in the 
figure) produces a 532 nm beam. The prism deflects the first harmonic (1064 nm) 
beam to a beam stopper. The 532 nm source is scattered by a rotating ground glass 
disk (the random media) in order to obtain the pseudo-thermal field. The beam 
is then split in three beams by means of two beam splitters (a polarizing beam 
splitter BS1 and a non-polarizing beam splitter BS2) and a 45° mirror and sent 
to a EMCCD camera. The presence of the half waveplate WP on the path of the 
beam before the rotating disk, together with the polarizing beam splitter, allows to 
balance the intensities on the three arms. 

The dimension of the coherence areas of the multimode thermal beams are set 
by the lens LI of focal length / = 750 mm. The lens is put in a / — / configuration 
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Fig. 2. Experimental set-up. 



respect to the pinhole and to the object in order to observe the far field speckle 
pattern at the plane with the object. In this way, a single transverse wavevector q is 
associated to a single point x = (A/ /2ir)q. The pinhole is at distance do = 110 mm 
from the rotating ground glass. As a matter of fact, the dimension of the coherence 
areas is also proportional, for a given wavelength, to the distance between the 
ground glass and the pinhole and inversely to the dimension of the pinhole. A 
second pinhole after LI determines the dimension of the source. 

The lens L2 before the camera is used to image the beams on the CCD. A 
magnification occurs in order to image the whole area of the three beams on the 
658 x 498 pixels sensitive area of the CCD, that is 6, 58 mm x 4, 98 mm since the 
pixel area is 10 /im x 10 /an. Distances di = 1420 mm and d 2 = 200 mm are chosen 
in order to obtain a proper demagnification. 

The camera is triggered with the pulses of the laser. The repetition rate of the 
laser is externally set to 500 Hz, while the exposure time of the camera is of 2 
ms. Since the pulse duration is no longer than 1 ns, each frame taken with the 
CCD contains only one pulse of the laser. The ratio between the repetition rate 
of the camera and the frequency of rotation of the disk avoids a periodic speckle 
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Fig. 3. The image of the three thermal beams from a single shot of the laser a) Correlated regions 
in arm 1 and 2. b)Correlated regions in arm 2 and 3. 

pattern detection. Moreover, the disk has to rotate slowly enough to be considered 
stationary with respect to the duration of the pulse of the laser. 
The object is an iron curled wire with a diameter of 1.5 mm. 



4. Toward third-order ghost imaging: experimental results 

Fig. ([3]) reports the image obtained from a single shot of the laser. We can appreciate 
the speckle structure and the correlations in its intensity on the three thermal 
beams. The object is in the central arm, named Si; S2 and S3 denotes the reference 
arms. In order to retrieve information about the correlations, the first step is to 
determine with high precision the regions R2 and R3 correlated to the one with the 
object, R\ (Fig. ([3]) a)). Due to the presence of the object, we first determine the 
distance between R\ and R2 using a region C\ without the object and finding 
the correlated one, C2 in the arm 2. Starting from a region 6*2,0 of S2 we displace 
the region in horizontal and vertical directions, pixel by pixel, in order to find the 
maximum of the normalized second order correlation coefficient pi for a frame k: 
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Fig. 4. a) The correlation coefficient as a function of the displacement vector D b) The peak of 
correlation. Its FWHM represents the diameter of the speckle 



c2 k = < JTC^TKg + d) > - < ii(xi) >< + d) > (6) 

(^l)V/^ (^2,0 + D) 

where l[ (x±) and 7^* (^2,0) ar e the intensities of the frame k for the region 
C\ and C2,o respectively. D is the displacement vector in pixel. In Eq. Q the 
symbol < .. > refers to spatial averages over the pixels of the selected regions. 
The distance between the correlated regions is computed as: di^ = \xi — with 
x% = x 2 fl + D max , where D max correspond to the displacement to be applied in 
order to reach the maximum of the correlation. Finally we determine the region i? 3 
correlated to R 2 using the same procedure (Fig. [3]b)). 

The correlation coefficient as a function of the displacement vector D is depicted 
in Fig. Q a) ; the part b) of the figure reports the section of the peak of correlation 
in the horizontal direction. The reported correlation between the regions is quite 
high (~ 0.88), denoting a low level of losses in the optical path and low noise. The 
FWHM of the peak is a good estimation of the size of a speckle; in our experiment 
we generate speckles of radius x co h ~ 3 pixels = 30 /Ltm. 

To perform ghost imaging, we use a bucket detector on arm 1. For each frame 
k taken with the camera, the measured quantity in arm 1 is then the integral over 
X\\ = J dxil[ (x\). In practice, it corresponds to the sum over the pixels of 
the region Ri with the object. The intensity distributions of the correlated regions 



14, 2010 0:19 WSPC/INSTRUCTION FILE MEDA2 



7 



r r t 




1 






" ft > 

i * 








Si 


a) 




c) 



Fig. 5. a) The image of the real object b) The second order ghost image c) The third order ghost 
image 



on the other arms are registered as an array of intensity values /• (xj ) . 

We reconstruct normalized second and third order correlation functions c2 and 
c3 as the ghost images of our real object. For c2 we used arm 2 and arm 3 alter- 
natively as reference. Fig. ^ shows our preliminary results. Part a) of the figure 
reports the single frame image of the object in the region Ri, while part b) shows 
one of the two reconstructed images by means of c2. Finally, the c) part is the 
third order ghost imaging of the object. The correlation coefficients are computed 
averaging over Nf = 400 frames. The reconstructed images are evident even if the 
resolution is quite low. This can be understood once you consider that the speckles 
have a diameter of 6 speckles, i.e. a dimension comparable to the one of the object 
(i.e. tens of speckles according to Fig. §5§ a)) 

As a figure of merit to estimate the improvement of third order ghost imaging 
respect to second order one we chose the visibility of the reconstructed images, here 
defined as: 

Y _ CJbgck — CJobj j>y\ 
CJback i CJXobj 

with j — 2, 3, where cjb ac k is the value of the correlation coefficient averaged over 
a region without the object and cj bj is computed averaging over a region with the 
reconstructed object. 

We obtain V2 = 0.2301 ± 0.0008 (the result is independent by the choice of arm 
2 or 3 as reference) and V3 = 0.269 ± 0.006. It is clear that there is an improvement 
in the value of the visibility of the reconstructed image by means of third order 
correlations. It was theoretically predicted by Bai et al. ^ that the existence of ad- 
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ditional correlation parts, namely correlations between the two reference detectors 
and between the three detectors, leads to an improvement of the visibility. Our 
experimentally measured increment in the value of V3 with respect to V2 can be 
compared to their theoretical predictions. Moreover, the gap remains stable even 
if we use different set of images. We also observe that the number of frames used 
to reconstruct the image is sufficient to reach the maximum gap. An increment of 
the number of frames do not lead to an improvement of the ghost imaging perfor- 
mances. In any case, we want to stress that the reported data are a preliminary 
experimental demonstration of the improvement in the visibility and that a more 
complete analysis is postponed to a subsequent paper. 

5. Conclusions 

In this paper we have presented a realization of third order ghost imaging of a real 
object exploiting a CCD camera. We also show an improvement in the visibility 
of third order ghost imaging respect to second order one. Having overcome some 
limits of previous experiments, it represents a relevant step toward real applications 
of ghost imaging E2I 
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